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Abstract. Stellar spectropolarimetry has become extremely popular during the last 
decade, and has led to major advances in our understanding of stellar magnetic fields 
and of their impact on stellar structure and evolution. Many important discoveries have 
been obtained thanks to observations performed with the FORS low-resolution spec- 
tropolarimeters of the ESO Very Large Telescope. We first review and summarise the 
major results of a homogeneous re-reduction and analysis of all single-slit FORS 1 spec¬ 
tropolarimetric observations. This work revealed a non-negligible dependence of the 
results upon the adopted reduction and analysis procedure, as well as the presence of 
instabilities, revealing that photon noise is not the only source of uncertainty. As a con¬ 
sequence of our new analysis and assessment of the uncertainties, we are not able to 
confirm a large number of magnetic field detections presented in the past for a variety 
of stars. We further summarise the results of FORS2 spectropolarimetric observations 
of the A0 supergiant HD 92207 which allowed us to explore further the nature of the 
instabilities, roughly constraining their maximum impact on the derived Stokes profiles 
and magnetic field values. We finally present new results obtained with a further inde¬ 
pendent pipeline on the FORS2 data of HD 92207, confirming our previous analysis, 
and discuss simple quality-check controls which can be performed on the data in order 
to distinguish between genuine and spurious signals. All together, our results reveal that 
the FORS spectropolari meters are indeed reliable instruments, when their capabilities 
are not pushed beyond the limits of a Cassegrain mounted low-resolution spectrograph. 


1. Introduction 

A large number of low resolution spectropolarimetric data aiming at detecting stellar 
magnetic fields have been collected with the FORS (FORS 1 and FORS2) spectropo- 
larimeters (Appenzeller et al. 1998) of the ESO Very Large Telescope (VLT). FORS1 
was decommissioned in April 2009 and the polarimetric optics were then moved to 
FORS2. Both FORS1 and FORS2 are multi-mode optical instruments capable of imag¬ 
ing, polarimetry, and long slit and multi-object spectroscopy. 
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Spectropolarimetric observations obtained with both instruments, particularly FORS1, 
have been used to claim the detection of magnetic fields in several different classes of 
stars across the Hertzsprung-Russell diagram, including for instance central stars of 
planetary nebulae (Jordan et al. 2005), hot subdwarfs (O’Toole et al. 2005), /jCcphci 
and slowly pulsating B stars (Hubrig et al. 2009a), Be stars (Hubrig et al. 2009b), and 
normal O-type stars (Hubrig et al. 2008, 2013). Many of these detections have been 
obtained close to the limit of the 3cr level and only some have been confirmed on the 
basis of further high resolution spectropolarimetric observations. 

As a matter of fact, in relation to FORS spectropolarimetric data there are incon¬ 
sistencies of various nature present in the literature: 

• inconsistencies between field measurements obtained with FORS and other in¬ 
struments (e.g., ESPaDOnS; see Silvester et al. 2009); 

• inconsistencies between the results obtained by different authors from the same 
FORS dataset (e.g., McSwain 2008; Hubrig et al. 2009b); 

• inconsistencies between the results obtained by the same author, but at different 
epochs, from the same FORS dataset (e.g., Wade et al. 2005, 2007); 

• inconsistencies between field measurements obtained from different subsets within 
a series of frames (e.g., Hubrig et al. 2004); 

• some global inconsistencies of the full FORS dataset (e.g., high incidence of field 
detections in the null profile; Bagnulo et al. 2012). 

In Sect. 2 we review the major results and conclusions obtained from a homoge¬ 
neous re-reduction and analysis of all single-slit FORS 1 spectropolarimetric observa¬ 
tions, first presented by Bagnulo et al. (2012). In Sect. 3 we further review the major 
results obtained from FORS2 spectropolarimetric data of the A0 supergiant HD 92207, 
presented by Bagnulo et al. (2013), disproving the detection of a magnetic field, as 
suggested by Hubrig et al. (2012). We then present the results of a further re-reduction 
of the FORS2 observations of HD 92207 using an independent pipeline (Sect. 3.1). In 
Sect. 3.2 we speculate upon the possible origin of a spurious Stokes V signal found in 
the FORS2 observations of HD 92207. To conclude, we discuss simple quality-check 
controls that can be performed on the data in order to distinguish between genuine and 
spurious signals, as well as suggestions aimed at obtaining the highest possible quality 
single-slit spectropolarimetric observations. 


2, The re-analysis of the ESO FORS1 spectropolarimetric data archive 


With the aim of quantifying these inconsistencies and in order to try to identify their 
origin we performed a re-analyses of all observations gathered in circular polarisation 
with the FORS 1 instrument. We performed the whole analysis using the same technique 
and tools, namely the ESO FORS pipeline (Izzo et al. 2010) for the data reduction 
and a suite of customised fortran codes for the measurement of the surface average 
longitudinal magnetic field «B,)) which was obtained using the following relation (see 
Angel & Landstreet 1970; Borra & Landstreet 1980): 


V(A) = -geffC z A 2 


1 d/(T) 


<B Z ) 


1(A) dA 
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and the least squares technique, first proposed by Bagnulo et al. (2002) and further 
refined by Bagnulo et al. (2012). In Eq. 1 V(A) and 1(A) are Stokes V and I, respectively, 
geff is the effective Lande factor, which was set to 1.25 except for the region of the 
hydrogen B aimer lines where g e ff was set to 1, and 


z 4nm e c 2 

where e is the electron charge, m e the electron mass, and c the speed of light (C, - 
4.67 x 10~ 13 A -1 G -1 ). Bagnulo et al. (2012) also give a discussion of the physical 
limitations of this technique. 

A typical polarimetric observations is composed of a sequence of spectra obtained 
alternatively rotating the quarter waveplate from -45° to +45° every second exposure 
(i.e., -45°, +45°, +45°, -45°, -45°, +45°, etc...). All details of the adopted data reduc¬ 
tion and analysis of the FORS1 archive are given in Bagnulo et al. (2012). We excluded 
from this re-analysis all data obtained in “multi-object” mode because of problems en¬ 
countered in the reduction of these data using the ESO FORS pipeline (see Bagnulo 
et al. 2012, for more details). 

The first results of the re-analysis are given in Bagnulo et al. (2012), while a further 
work (Bagnulo et al. in prep) will present the magnetic field values we obtained from 
the whole FORS1 dataset. In Bagnulo et al. (2012) we gave a re-assessment of the 
most controversial magnetic field detections claimed in the past on the basis of FORS 1 
spectropolarimetric data. In this work we concluded also that the data reduction and 
analysis procedure has a large impact on both magnetic field and uncertainty values and 
that the data are occasionally affected by systematics which arc not taken into account 
in the error budget. 

To present the effects of adopting different reduction procedures on the finally ob¬ 
tained (B ; ) values we showed the results gathered for three stars: a clearly magnetic 
chemically peculiar A-type star (HD94660), a non-magnetic star (HD96441), and a 
star consistently presenting a magnetic field detection around the 3 cr level (HD 171184), 
where for each star we applied twelve slightly different reductions and analysis proce¬ 
dures (see Table 2 of Bagnulo et al. 2012). The results we obtained for the latter star 
(i.e., HD 171184) are of particular' interest: out of the twelve (B z ) values we derived, 
four show a magnetic field detection at less than 3cr, seven show a magnetic field de¬ 
tection between 3 cr and 4 cr, and one shows a magnetic field detection larger than 4 cr. 
Here we show that 3-to-4cr magnetic field detections could be turned into less than 3<x 
detections simply by e.g., trimming the first and last 10 data-points in Stokes I and V, 
or by applying a 2 pixel rebinning (see Bagnulo et al. 2012, for more details). 

In this respect, the major problem is that it is not possible to clearly identify the 
“best”/“most appropriate” data reduction and analysis procedure. For example, one 
could decide to do a background subtraction or not and it is not clear which of the 
two options is the correct one, but they would lead to two different (B-) values, where 
the difference is not quantifiable a priori ; the situation is similar, for example, for flat 
fielding. One could also decide to apply a sigma clipping or not, and there are dif¬ 
ferent ways of doing it, all of them equally valid and probably leading to (possibly 
significantly) different (B ; ) values. 

In Bagnulo et al. (2012) we further performed a statistical analysis of the magnetic 
field values obtained from the diagnostic null profiles (N Donati et al. 1997; Bagnulo 
et al. 2009), hereafter denoted as (N,), which give and indication of the noise level and 
are used to highlight the presence/absence of systematic errors in the data (see Bagnulo 
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et al. 2009, for more details). In Fig. 7 of Bagnulo et al. (2012) we presented the 
histogram of the (B-) values obtained from the null profiles in the whole FORS1 catalog 
normalised to their uncertainties. The histogram should take the form of a Gaussian 
with cr=l, but we obtained something which resembles a Gaussian with cr= 1.365; this 
clearly indicates the presence of systematics in the FORS 1 data and that photon noise is 
not the only source of uncertainty. This implicitly explains, at least partly, the number 
and nature of inconsistencies listed in Sect. 1 

These tests and considerations led us to the conclusion that the 3cr level cannot be 
considered as an appropriate detection level for magnetic field measurements obtained 
with the FORS spectropolarimeters, though this consideration extends in general to 
all long-slit low resolution spectropolarimeters. On the basis of our large re-analysis 
and by comparing results obtained with different instruments we concluded that a solid 
magnetic field detection with the FORS spectropolarimeters can be established only 
in the presence of repeated 5-to-6cr detections: as shown below, even a single clear 
detection (>5cr) might not be sufficient to safely avoid spurious detections. 


3. FORS2 spectropolarimetric observations of HD 92207 

Hubrig et al. (2012) presented the results gathered from the analysis of FORS2 spec¬ 
tropolarimetric observations of the A0 supergiant HD 92207 obtained in three different 
epochs with the grism 600B. Their analysis showed the detection of a magnetic field at 
two different epochs: the first one at the ~9cr level and the second one at the ~4cr level. 
On the basis of the conclusions given by Bagnulo et al. (2012) this detection might ap¬ 
peal - to be genuine. Nevertheless, the Stokes V profiles shown by Hubrig et al. (2012) 
revealed the presence of several suspiciously strong Stokes V signals at the position of 
sharp and deep lines (i.e., hydrogen Balmer lines and CanH&K lines) which appeal - to 
be too large for the detected magnetic field strength. For this reason we retrieved the 
FORS2 data from the ESO archive and re-analysed them in order to confirm/disprove 
the magnetic field. 

3.1. The spurious detection of a non-existing magnetic field 

We performed the data reduction and analysis using two independent tools. The first 
one is that used for the re-analysis of the FORS1 archive (see Sect. 2 and Bagnulo et al. 
2012), while the second one is based on a suite of custom made iraf 1 (Tody 1993) and 
idl routines which follow in part the technique and recipes presented by Bagnulo et al. 
(2012, 2013). 

Here we briefly describe the major characteristics of the latter reduction and anal¬ 
ysis procedure/tools while more information will be given in a separate paper (Fossati 
et al., in prep.). The iraf pipeline reduces the raw frames only by removing the bias, i.e. 
without correcting for the flat-field, and it performs an average extraction of the spectra 
without background subtraction. Within each dataset, each parallel/perpendicular beam 
is wavelength calibrated using the parallel/perpendicular beam of one wavelength cal¬ 
ibration lamp, usually that obtained with the position angle at -45°. To make sure that 


'Image Reduction and Analysis Facility (IRAF - http://iraf.noao.edu/) is distributed by the Na¬ 
tional Optical Astronomy Observatory, which is operated by the Association of Universities for Research 
in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 
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different sets of arc lines arc not used for the wavelength calibration, the iraf pipeline 
requests human interaction to calibrate the parallel and perpendicular beams in a com¬ 
pletely independent way and it makes sure that always the same arc lines and the same 
fitting function (i.e., a 3 rd order spline) is used for the two beams. After having ap¬ 
plied the wavelength calibration, the spectra arc binned with the natural sampling of 
the instrument, 0.75 A/pix in the case of the FORS2 data of HD 92207. 

The extracted and wavelength calibrated spectra arc combined by an idl routine 
in order to obtain Stokes /, V, and the diagnostic N parameter using the “difference” 
method following the formalisms of Bagnulo et al. (2009). One can optionally use the 
uncertainties given by the iraf extraction package (apall) or opt for pure photon noise 
error bars; the latter is the default setting as we noticed that the default iraf uncertainties 
tend to be systematically smaller. The Stokes V profile can be optionally rectified; here 
we used a 4 th order polynomial. A sigma clipping can also be applied on the basis of 
the N parameter; here we filtered out all data-points where the N profile deviates more 
than 3<x from N, where cr is the standard deviation of N. Both (B.) and (N.) values arc 
finally calculated by minimising the expression given in Eq. 7 of Bagnulo et al. (2012) 
using either the hydrogen lines, or the metallic lines, or the whole spectrum within the 
observed wavelength region. The (B-) and (N,) uncertainties arc then rescaled by the 
X 2 as suggested by Bagnulo et al. (2012). 

Bagnulo et al. (2013) presented the results of the re-analysis of the FORS2 spectra 
of HD 92207 that showed the largest (B-) value in the analysis conducted by Hubrig 
et al. (2012) (i.e., data obtained in 2011). Using the reduction and analysis procedure 
shown in Bagnulo et al. (2012) they obtained a 2.9cr magnetic field detection from the 
Stokes V profile and a 5.7 cr detection from the N profile. This indicates that a magnetic 
field signature is not present in the analysed FORS2 data, in contrast to the results of 
Hubrig et al. (2012). In particular, the clear (N,) detection indicates the presence of 
systematics in the data, the origin of which will be discussed in Sect. 3.2. 

As a further check of this result, we re-analysed the 2011 FORS2 data using the 
iraf/idl reduction and analysis tools described above. Figure 1 shows the graphical out¬ 
put of the iraf/idl pipeline obtained when considering the hydrogen lines. The results 
are given Table 1. These independent results arc in agreement with those obtained by 
Bagnulo et al. (2013) and in clear disagreement with those obtained by Hubrig et al. 
(2012). For the analysis performed with the iraf/idl we used the same wavelength range 
shown by Hubrig et al. (2012) and Bagnulo et al. (2013). 


Table 1. (B-) and (N,) values obtained from the FORS2 data of HD 92207 using 

the iraf/idl pipeline and the while spectrum (2nd and 3rd columns) or the hydrogen 
lines (4th and 5th columns). 


MJD 

(B z )all 

[G] 

(Nj)all 

[G] 

(Bj)hyd 

[G] 

-o 

m i — i 

"bO 

55688.168 

-249+83 

-281+78 

-182+99 

-254+93 

55936.341 

25+129 

-210+119 

18+143 

-232+140 

56018.224 

-68+158 

79+169 

95+182 

8+192 
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These results show also that the more “manual” iraf/idl reduction leads to smaller 
(N,) detections compared to the analysis computed with the fully automatic ESO FORS 
pipeline. Such a trend is also present for a few other FORS 1 data which showed rather 
large (N,) detections in Bagnulo et al. (2012) and that we re-analysed with the iraf/idl 
pipeline as double-check. We believe that the difference lies in the way the wavelength 
calibration is performed within the two pipelines. The manual wavelength calibration 
performed with the iraf/idl pipeline allows one to make sure the same set of lines is 
used for both parallel/perpendicular beams, in spite of the fact that this leads to the use 
of a limited number of calibration lines (usually 6-7 for each beam with the grisms 
600B, 600V, and 1200B). On the other hand, the fully automatic ESO FORS pipeline 
aims instead at maximising the number of lines used for the wavelength calibration, but 
this often leads to the use of a slightly different set of calibration lines for the two beams, 
because of their slight difference in transmission/efficiency (i.e., the weakest lines used 
for the calibration of one beam become too weak to be used for the calibration of the 
other beam). We are still exploring whether this effect might actually be at the origin of 
the anomalously large number of (N z ) detections in the FORS1 archive (Bagnulo et al. 
2012 ). 

We also re-analysed with the iraf/idl pipeline the FORS2 data of HD 92207 ob¬ 
tained in 2012, consistently obtaining non-detections, as shown by the results given in 
Table 1. We can therefore firmly conclude that there is no evidence of the presence of 
a magnetic field in the FORS2 data of HD 92207. 

Bagnulo et al. (2013) also presented the results of a magnetic field search for 
HD 92207 conducted with the HARPspol spectropolarimeter (Snik et al. 2011; Piskunov 
et al. 2011) feeding the HARPS spectrograph (Mayor et al. 2003) attached to the ESO 
3.6-m telescope in FaSilla, Chile. The data, obtained in 2013 on four different epochs, 
were reduced using the reduce package (Piskunov & Valenti 2002) and analysed using 
the Feast-Squares Deconvolution technique (FSD; Donati et al. 1997; Kochukhov et al. 
2010). The analysis of the FSD profiles led to clear non-detections with uncertainties on 
(B-) of the order of 10 G. It is important to notice that, because of the cross-over effect, 
the presence of a structured magnetic field, particularly of the strength suggested by the 
results of Hubrig et al. (2012), would be clearly visible in the FSD profiles regardless 
of the rotation phase. These results firmly indicate that there is no trace of a strong 
(>100 G) large-scale magnetic field in HD 92207, in contrast to what was concluded by 
Hubrig et al. (2012). 

Despite our clear results on the non-detection of a magnetic field for this star, it 
is important to explore the origin of the spurious detection obtained by Hubrig et al. 
(2012): as a matter of fact, they obtained a ~9cr magnetic field detection from the 2011 
FORS2 data, well above the “safety” threshold of 5-6 cr suggested by Bagnulo et al. 
(2012). We carefully inspected all major steps performed during the data reduction in 
order to identify the possible origin of the spurious signal obtained by Hubrig et al. 
(2012). We were finally able to reproduce Hubrig et al. (2012)’s results by deliberately 
doing an inaccurate wavelength calibration (see Bagnulo et al. 2013, for more details). 
It is important to remark that by deliberately doing an inaccurate wavelength calibration 
we simply introduced a systematic wavelength shift in the parallel beam relative to 
the perpendicular beam. The Stokes V profiles we obtained in this way (see Fig. 6 
of Bagnulo et al. 2013) present spikes at the position of narrow and deep lines, very 
similar to those shown by Hubrig et al. (2012). From the analysis of the Stokes V 
profile we obtained (B-) = -325±105 G, again very similar to that obtained by Hubrig 
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et al. (2012), except for the uncertainty. Our uncertainty was larger because we rescaled 
the error bar with the y 2 °f d lc linear fit. We double-checked the procedure and results 
described here using the iraf/idl, obtaining the same profiles and arriving at the same 
conclusions as in Bagnulo et al. (2013), as shown in Fig. 2 and Fig. 3. 

Nevertheless, in contrast to what obtained by Hubrig et al. (2012), we obtained 
similar spikes in the N profile as well, which led to (N-) = -355±75 G, clearly indicat¬ 
ing the spurious nature of the spikes. This result is somewhat reassuring as it shows 
that the N profile has “done its job” by clearly indicating that the Stokes V detection 
should not be considered to be reliable. Similarly, by rescaling the uncertainty with the 
X 2 of the lineal - fit the detection level obtained from both V and N decreased consider¬ 
ably, showing that the “safety” threshold of 5-6 cr suggested by Bagnulo et al. (2012) 
is indeed valid. 

We finally attempted to reproduce also the rather flat IV profile shown by Hubrig 
et al. (2012). In the end, by reshuffling the order of the frames used to derive the 
Stokes parameters and the null profile, we could obtain a much smoother N spectrum, 
qualitatively si mi lar to that presented by Hubrig et al. (2012). In other words, given a 
set of spectropolarimetric observations, the N profile is not unique (see Bagnulo et al. 
2013, for more details). 

3.2. On the origin of the short-time line shifts observed in the FORS2 spectra of 
HD 92207 

In Bagnulo et al. (2013) we used the CauK line to show the presence of line shifts, 
of the order of about a quarter of a pixel (~4 km s -1 ), in spectra obtained in consecu¬ 
tive frames. We also noticed that such shifts are consistently present across the whole 
spectrum, but concentrated on the analysis of the CauK line because of its strength 
and depth. We considered a range of possible astrophysical and systematic origin of 
these shifts, but finally excluded all astrophysical (e.g., pulsation) nature of the shifts 
and attributed them to systematics caused by either instrument flexures or observing 
conditions (e.g., seeing variations). We showed that these systematic shifts are at the 
origin of /) the slight spurious signal which is consistently present in both (B-) and (N-) 
values obtained with both our data reductions and it ii) a clear spike in the Stokes V 
profile at the position of the CauK line (see Bagnulo et al. 2013, for more details). 

Hubrig et al. (2014) confirmed what was found by Bagnulo et al. (2013) in terms 
of line shifts, but attributed them to physical changes in the star suggesting short-time 
pulsations as the most likely cause. To support this conclusion they showed the lack 
of similar line shifts for two other stars (HD 93843 and ( Oph) observed in the same 
nights in which the FORS2 spectra of HD 92207 were obtained. 

It is here important to recall that each FORS2 observation of HD 92207 at each 
position angle was performed with an exposure time of 3 seconds and that between 
two consecutive exposures in a sequence there is about 1 minute (CCD readout). This 
implies that, by assuming the line profile variations detected by Bagnulo et al. (2013), 
and then confirmed by Hubrig et al. (2014), find their origin in the stellar pulsation, the 
period has to be of the order of a few minutes, at most. Pulsating supergiants, similar 
to HD 92207, have typical minimum pulsation periods of the order of days, while for 
giants, dwarfs, and white dwarfs the typical minimum pulsation period of is of the order 
to hours, minutes, and seconds, respectively. In other words, it is extremely unlikely, 
if not even physically impossible, that an A0 supergiant such as HD 92207 presents 
pulsations with a period typical of that of main-sequence stars or white dwarfs. In 
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addition, it is really hard to imagine how a star, with a radius of the order of 140 R Q 
(Hubrig et al. 2012) can show pulsation with a period of the order of minutes; note that 
a free photon moving on a straight line in the void would need about 5.4 minutes to 
cover a distance of 140 R Q . 

The presence of similar line shifts for all observed lines and the physical character¬ 
istics of the star (e.g., radius, rotation rate, temperature, lack of a strong magnetic field) 
tends to exclude further astrophysical origins, such as surface spots, wind clumping, 
etc ..., of the observed line shifts. As a result, we have to search for an instrumental or 
observational systematic origin of the line shifts. 

Let us examine the exposure times, adopted slit widths, seeing conditions, and 
coherence times 2 of the FORS2 observations of HD 92207, HD 92843, and £ Oph, the 
main subject of this work and the two stars used by Hubrig et al. (2014) to support their 
claim that the observed line shifts arc of astrophysical origin. Table 2 lists the relative 
quantities of each star. 


Table 2. Columns one and three give respectively the slit width (in ") and 
exposure times (in seconds) adopted for the FORS2 observations of HD 92207, 
HD 92843, and Oph. Columns two and four give the sky conditions, namely see¬ 
ing (in ") and coherence time (in seconds) at the time of the FORS2 observations of 
HD 92207, HD 92843, and f Oph and listed in the fits headers. 


Star 

name 

Slit 

width 

["] 

Seeing 

["] 

Exposure 

time 

[s] 

Coherence 

time 

[s] 

HD 92207 

0.4 

~0.5 

3 

6 

HD 92843 

0.4 

-0.5 

30 

6 

<TOph 

0.4 

-0.7 

0.25 

4 


For HD 92207, the FORS2 observations have been performed using a narrow slit 
and in excellent seeing conditions, namely with a seeing similar to the slit width. In 
addition, the adopted exposure time was of the same order of the timescale of seeing 
variations. The most likely explanation of the observed line shifts lies therefore in vari¬ 
ations of the position of the star on the slit caused by seeing variations in concurrence 
with the adoption of a slit width similar to the seeing conditions. From the values 
listed in Table 2 it appeal's then clear why similar line shifts are not observed for both 
HD 92843 and ( Oph: for the former the adopted exposure time is much greater than 
the timescale of seeing variation, while for the latter the adopted exposure time is much 
shorter than the timescale of seeing variations. For the observations of both HD 92843 
and £ Oph the fact that the adopted slit width is comparable to the seeing does not play 
a role because of the different timescales between the exposure times and the coherence 
times. The observed line shifts are most likely also at the origin of the spurious (N z ) 
detections reported in Sect. 3.1. 


2 The coherence time is given in the fits header and gives an indication of the timescale of seeing variations. 
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4 . Conclusions 


We re-reduced and analysed the whole FORS1 ESO archive of long-slit spectropolari¬ 
metric observations concluding that a magnetic field detection can be safely reported 
only on the basis of repeated 5-6 cr detections. In addition, we reported the presence of 
excess noise in the FORS1 data, that we ascribe to systematics of multiple origin (e.g., 
instrumental, observing conditions, etc ...). We also realised that different, but equally 
valid, data reduction procedures and tools lead to (significantly) different results (see 
Bagnulo et al. 2012, for more details). 

We further proceeded to an extensive and deep analysis of the FORS2 observa¬ 
tions of HD 92207 that led Hubrig et al. (2012) to report the detection of a rather strong 
magnetic field ((B,) —400 G) at the ~9cr level. We re-analysed the FORS2 data us¬ 
ing two independent pipelines consistently obtaining non-detections. We managed to 
reproduce the strong spikes clearly visible in the Stokes V spectra shown by Hubrig 
et al. (2012) only by deliberately doing an inaccurate wavelength calibration that intro¬ 
duced a systematic shift between the parallel and perpendicular beams. Despite this, 
we obtained a null profile that clearly indicated the spurious nature of the Stokes V 
spikes. In contrast, Hubrig et al. (2012) showed smooth N spectra, that we were able to 
reproduce by reshuffling the order of the frames used to derive the Stokes parameters 
and the null profile (see Bagnulo et al. 2013, for more details). We reported also a non¬ 
detection of a magnetic field at the level of a few tens of gauss using high resolution 
spectropolarimetric HARPSpol data analysed with the FSD technique. 

In Bagnulo et al. (2013) we also showed the presence of systematic shifts in the 
line profiles of Stokes I spectra obtained from consecutive frames. The observed line 
shifts have been further confirmed by Hubrig et al. (2014) who attributed them to stellar 
pulsations. If pulsation is really present and if it causes the observed line profile varia¬ 
tions, the pulsation period has to be of the order of a few minutes (5 minutes at most). 
Such a short pulsation period is typical of dwarfs and white dwarfs and certainly not of 
supergiants with stellar radii larger than 100 R Q . Further characteristics of the star, such 
as the lack of a strong magnetic field, the rotation period, and the effective tempera¬ 
ture, exclude other astrophysical origin for the line profile variations. In Bagnulo et al. 
(2013) we concluded then that one has to seek their origin in the instrument flexures or 
in the observing conditions. 

As a matter of fact, the observations of HD 92207 have been obtained with an 
exposure time comparable to the timescale of seeing variations. This, in addition to the 
adoption of a slit width similar to the seeing, has led to a non-negligible instability of 
the star’s image on the slit and therefore to the observed line shifts. 

We have to remark that the systematics observed for HD 92207 appeal - to be of 
a rather extreme magnitude. Nevertheless, by following our reduction procedures we 
systematically obtained both a field detection smaller than what recommended by Bag¬ 
nulo et al. (2012) for a reliable detection and a significant signal in the null profile, 
indicative of the presence of systematics. 

On the basis of our experience we recommend to always re-scale the uncertainties 
of the (B-) and (N,) measurements by the y 2 value and to report a reliable magnetic 
field detection only in presence of a repeated 5-6 cr detection (Bagnulo et al. 2012). We 
recommend to always look for the absence of line profile shifts in spectra obtained at the 
same position angle within a single sequence of polarimetric observations similarly to 
what done by Bagnulo et al. (2013). We also suggest to check the shape of the Stokes V 
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profiles obtained from the single couples of -457+45° frames for various lines, paying 
particular attention to the narrow and deep lines. This will clearly and immediately 
reveal the presence of systematics which will certainly affect the (B,) measurements. 
We finally recommend to always avoid very short exposure times, particularly if the 
observations arc performed with a slit width of comparable size (or larger) than the 
seeing conditions. 

To conclude, our work shows that, assuming one follows a basic set of observa¬ 
tional and analysis recommendations, it is indeed possible to confidently use instru¬ 
ments, such as FORS, to measure stellar magnetic fields, even of modest strength. 
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Figure 1. Graphical output of the iraf/idl pipeline obtained from the analysis 
of the 2011 FORS2 data of HD 92207 considering the hydrogen lines. Top panel: 
derivative of Stokes I. The regions used for the calculation of the magnetic field are 
marked by a thick blue line centered at zero. Middle panel: the top profile shows 
Stokes I arbitrarily normalised to the highest value, the middle red profile shows 
Stokes V (in %), while the bottom blue profile shows the N profile (in %). The 
green asterisks present on the Stokes V and N spectra mark the points which have 
been removed by the sigma clipping. The pale blue strip drawn behind the N profile 
shows the uncertainty associated with each spectral point. As a double-check of the 
statistical consistency of the data, the thick green bar on the left side of the N spec¬ 
trum shows the standard deviation of the N profile. Bottom-left panel: linear fit used 
for the determination of the magnetic field value using Stokes V (i.e., (B,)). The 
red solid line shows the best fit. From the linear fit we obtained (B,) = -182+99 G. 
Bottom-right panel: same as the bottom-left panel, but for the null profile (i.e., (N-)). 
From the linear fit we obtained (NT) = -254+93 G. 
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Figure 2. Same as in Fig. 1 but deliberately doing an inaccurate wavelength cal¬ 
ibration as described by Bagnulo et al. (2013). From the linear fits we obtained 
<B ; > = -242+211G and (N-) = -173+121G. 
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Figure 3. Same as in Fig. 2 but reshuffling the order of the frames used to de¬ 
rive the Stokes parameters and the N profile as described by Bagnulo et al. (2013). 
From the linear fits we obtained (B ; ) = -440+247 G and (N-) = -158+89 G. Note 
that by not rescaling the uncertainties on (B,) and (N,) by the v 1 value we obtained 
<B-> = -440+74 G and (N-> = -158+74 G. 











































